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Solar  Wind  Study  by  the  Kinematic  Simulation  Method 

During  the  last  12  months,  we  have  developed  a  3-D  kinematic 
simulation  method  in  order  to  be  able  to  deal  realistically  with  solar 
events.  In  particular,  we  have  developed  the  computer  code  in  order  to 
be  able  to  handle  the  propagation'  of  shock  waves  generated  by  six 
successive  solar  flares  from  the  same  active  region.  We  constructed 
first  a  steady  state  corotating  pattern  for  the  so-called  ‘two-sector’  or 
1  two-stream*  situation.  For  this  purpose  the  distribution  of  the  solar 
wind  speed  on  the  source  surface  (a  spherical  surface  of  radius  of  2.5 
Re)  is  assumed  to  have  the  minimum  speed  (V  ■  300  km/sec)  along  the 
heliomagnetic  equator,  which  is  assumed  to  be  Inclined  by  20°  with 
respect  to  the  heiiographic  equator;  the  speed  is  assumed  to  increase 
toward  higher  latitudes  in  both  the  northern  and  southern  hemispheres. 
In  this  way,  a  fast  and  a  slow  solar  wind  flow  stream  out  alternatively 
along  a  radial  line  fixed  in  space. 

The  first  step  of  our  method  is  to  construct  an  initial  distance  (R) 
-  time  (t)  relationship  for  particles  streaming  out  along  a  (solar) 
radial  line  fixed  in  space.  From  two  R-t  curves  at  two  successive 
instants  t^  and  t2  (t2  -  t^  *  At),  we  construct  the  velocity  (V)  - 

distance  (R)  relationship.  Although  our  method  is  basically  kinematic, 
faster  particles  are  not  allowed  to  take  over  slower  particles.  This 
restriction  alters  the  initial  R-t  relationship.  We  modify  the  initial 
R-t  relationship  until  the  resulting  V-R  relationship  resembles  closely 
an  empirically  constructed  V-R  relationship  on  the  basis  of  space  probe 
observations  or  of  theoretical  results;  for  the  R-t  and  V-R 


relationships,  see  Figures  i.4a,  1.4b  and  1.5  in  Hakamada  and  Akasofu 
(1982).  The  final  R-t  relationship  thus  obtained  becomes  the  basis  for 


our  study.  The  magnetic  field  Is  compressed  In  the  shocks;  the  flux 
density  and  field  Intensity  can  be  computed  using  flux  conservation.  The 
resulting  magnetic  field  configuration  in  the  equatorial  plane  is  the 
familiar  corotating  structure  which  consists  of  two  'spiral  arms'  (the 
corotating  interaction  region). 

Effects  of  solar  flares  are  introduced  into  this  steady  corotating 
pattern  by  adding  a  high  speed  flow  from  a  circular  area,  centered  around 
a  solar  flare.  In  the  circular  area,  the  flow  speed  is  assumed  to  have  a 
Gaussian  distribution  of  the  'half-width'  a;  the  flow  speed  Vp  at  the 
center  is  assumed  to  vary  with  time  as  Vp  »  ^Fmax  *  *  e  thus,  a 

flare  is  characterized  by  six  parameters,  the  latitude  (9),  longitude  (4>) 
and  the  onset  time  (Tp)  of  a  flare,  the  maximum  flow  Vpmax  (km/sec)  at 
the  center  of  the  circular  area  of  size  parameterized  by  a(°),  and  the 
time  variations  of  Vpfflax  by  x( hr). 

The  first  flare  is  assumed  to  occur  at  T  >  0  and  at  0  *  0° ,  <t>  * 
0®.  This  particular  flare  is  parameterized  by  the  maximum  wind  speed 
Vpmax  *  800  km/sec,  x  -  12  hrs  and  o  =*  60®.  The  second  flare  is 

introduced  48  hours  (2  days)  after  the  first  flare.  It  is  assumed  that 
the  same  active  region  is  responsible  for  the  second  flare;  at  the  time 
of  the  second  flare,  the  active  region  rotated  by  =  28.3®  from  the  time 
of  the  first  flare.  The  second  flare  is  assumed  to  be  characterized  by  a 
faster  flow  Vpmax  *  1000  km/sec,  the  same  growth-decay  curve  (x  *  12  hrs) 
and  a  narrower  extent  (o  -  10°)  than  the  first  one,  and  so  on. 

Figure  1  shows  the  disturbance  patterns  in  the  ecliptic  plane  at 
distances  <  5  au  at  T  ■  4.0,  5.0,  6.0  and  7.0  days  after  the  onset  of  the 
first  flare.  At  T  ■  4.0  days,  the  one  can  see  that  the  familiar  spiral 
structure  of  the  interplanetary  magnetic  field  with  the  two  'spiral 
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region.  The  spiral  curves  Indicate  the  magnetic  field  lines  In  the  ecliptic  plane.  The  figure  shows  the 
propagating  structure  at  4,0,  5.0,  6.0  and  7.0  days  after  the  onset  of  the  first  flare.  The  shock  waves 
cause  narrow  belts  of  concentrated  magnetic  fluxes.  The  outer  limit  of  each  circular  area  Is  at  5  au. 


arms'  (  but  one  of  them  Is  disturbed  by  the  two  shock  waves.  The  third 
flare  began  at  T  =*  4.0  days,  but  the  associated  shock  wave  is  too  close 
to  the  sun  at  that  time. 

At  T  ■  5.0  days,  the  shock  wave  caused  by  the  third  flare  becomes 
evident.  The  fourth  flare  occurred  at  that  time.  At  T  =*  6.0  days,  the 
shock  wave  associated  with  the  fourth  flare  becomes  evident.  The  second 
shock  wave  catches  up  with  the  first  shock  wave  because  the  former  is 
associated  with  a  faster  flow  than  the  second  flow. 

At  T  *  7.0  days,  the  first  and  second  shock  waves  reach  one  of  the 
spiral  arms  at  a  distance  of  ~  4  au.  Note  that  the  two  spiral  arms 
rotate  steadily  with  the  sun  during  the  entire  period.  Note  also  that 
the  fifth  flare  begins  at  that  time.  Figure  2  shows  the  disturbance 
pattern  within  a  radial  distance  of  <  15  au,  at  T  *  7,  10,  15,  17,  19  and 
25  days  after  the  onset  of  the  first  flare.  At  T  *  10  days,  the  last 
(sixth)  flare  took  place. 

One  can  easily  see  that  the  shock  waves  and  the  corotating  spiral 
arms  interact,  forming  a  large-scale  magnetic  structure  and  expanding 
rapidly  outward.  The  corotating  structure  in  the  inner  heliosphere  was 
destroyed  by  the  successive  shock  waves.  However,  after  the  last  flare 
which  took  place  at  T  -  10  days,  it  begins  to  reform  near  the  sun.  It  is 
clearly  seen  at  T  *  25  days;  see  the  arrow  in  the  figure. 

McDonald  et  al.  (1981)  observed  a  large  decrease  of  galactic  cosmic 
ray  intensity  (in  the  outer  heliosphere)  a  few  months  after  successive 
flare  activities  in  April  1978  and  March,  April  and  May  1979.  As 
suggested  by  Burlaga  (1982)  and  Intrlllgator  and  Miller  (1982),  it  is 
quite  likely  that  the  flare-induced  shocks  and  corotating  structure  have 
profound  effects  on  the  propagation  of  both  galactic  and  solar  cosmic  ray 


15  AU  15  AU 


Figure  2.  Solar  wind  disturbances  caused  by  six  successive  flares  during  a 
period  of  10  days.  The  figure  shows  the  propagating  structure  at  T  *  7.0, 
10.0,  17.0,  15.0,  19.0  and  25.0  days  after  the  onset  of  the  first  flare.  The 
outer  limit  of  each  circular  area  is  at  15  au. 


particles  in  the  outer  heliosphere  by  forming  a  sort  of  barrier  to 
them. 

It  should  be  noted  that  the  present  study  is  a  kinematic  one  in 
which  a  simple,  empirical  method  (Hakamada  and  Akasofu,  1982)  is  used  to 
provide  a  first  order  construction,  temporally  and  spatially,  of  flare¬ 
generated  shocks  and  their  multiple  interactions  with  each  other  as  well 
as  with  corotating  interaction  regions.  It  is  hoped  that  such  a  first 
order  effort  will  be  of  some  use  in  interpreting  solar  wind  and  cosmic 
ray  observations  by  space  probes.  Ue  have  not  simulated  the  other 
dynamic,  thermodynamic,  and  magnetic  properties  (other  than  first-order 
IMF  distortion)  that  can  be  found  only  from  the  MHD  solutions  such  as 
those  by  Wu  et  al.  (1979,  1983),  Dryer  et  al.  (1980)  and  D'Uston  et  al. 
(1981). 

Figure  3  shows  the  computed  V-t  relationship  for  the  situation 
described  in  the  above.  It  may  be  compared  with  an  observation  of  the 
solar  wind  by  two  space  probes,  Voyager  1  (~  8  au)  and  Helios  1  (~  0.5 
au)  in  Figure  4.  The  velocity  profiles  at  1.0  and  10.0  au  in  Figure  3 
are  similar  to  those  observed  by  the  Helios  and  Voyager,  respectively. 
These  and  other  results  will  be  published  in  the  following  two  papers. 

S.-I.  Akasofu  and  K.  Hakamada,  Solar  wind  disturbances  in  the  outer 
heliosphere,  caused  by  six  successive  solar  flares  from  the  same 
active  region,  Geophys.  Res.  Lett,,  (in  press),  1983. 

S.-I.  Akasofu,  K.  Hakamada,  and  C.  Fry,  Solar  wind  disturbances  caused 
by  solar  flares:  Equatorial  plane,  Planet.  Space  Sci.,  (in  press), 
1983. 
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Figure  1.  The  time  variations  of  the  solar  wind  and  the  IMF  magnitude  at 
radial  distance  of  1,  5  and  10  au  along  the  radial  line  of  longitude  “  45 
•Vo  is  reckoned  counterclockwise  from  the  righthand  horizontal  axis. 


Figure  4.  The  solar  wind  speed  observation  by  Voyager  l  at  a  radial  distance  of  7.5  ~  8.5  au  and  by 
Hellos  l  at  a  radial  distance  of  +  0.5  au  (with  a  time  delay  assuming  a  constant  speed  of  500  km/sec);  after 

Burlaga  et  al.  (1983). 


At  present  we  are  studying  the  following  two  events: 


1 .  November  22  -  December  6,  1977 

This  was  a  unique  period  during  which  data  from  six  spacecraft 
are  available,  and  the  observation  was  summarized  in  detail  by 
Burlaga  et  al.  (1979)  in  their  paper  entitled  Interplanetary 
particles  and  fields,  November  22  -  December  6,  1977:  Helios, 
Voyager,  and  Imp  observations  between  0.6  au  and  1.6  au." 

2.  June  -  July  1,  1982 

There  were  three  active  regions  which  produced  a  large  number 
of  flares.  There  are  deep  space  probe  data  (Pioneer  10  and  11) 
during  this  period. 

II.  Calibration  of  the  Kinematic  Method  of  Studying  the  Solar  Wind  on  the 
Basis  of  a  One-dimensional  MHD  Solution 

The  empirical  kinematic  method  (the  HA  method)  can  be  Improved  by 
having  a  theoretical  or  observational  velocity  (V)  -  distance  (R) 

relationship.  Thus,  it  is  Important  to  make  every  effort  to  improve  the 
HA  method  on  the  basis  of  MHD  solutions  (the  V-R  relationship).  Once  the 
improvement  is  made,  the  HA  method  becomes  a  powerful  tool  in  inferring 
the  solar  wind  flow  and  the  interplanetary  magnetic  field  conf iguration 
in  the  heliosphere.  We  shall  show  some  of  the  new  results  which  are 
obtained  using  the  new  V-R  relationship. 

We  assume  that  the  interplanetary  medium  can  be  represented  a3  a 
single  fluid  with  negligible  dissipation.  In  terms  of  the  spherical 
coordinates  (r,  9,  $)  the  dependent  variables  are  density  p,  thermal 
pressure  p,  fluid  velocity  (Vr,  0,  V^).  The  independent  variables  are 


\ 

J  . 
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radius  r  and  time  t.  With  these  assumptions,  and  using  MKS  units  the  MHD 
equations  to  be  solved,  In  the  ecliptic  plane,  are  as  follows: 


ST  (f2p)  +  S7  <r2  PV  *  0 


a  ,  2  „  .  ^  a  2  „  2  .  .  <j> 

at  (r  PV  +F7  r  pVr  +  p  +  2? 


2rp  +  rpV^  -  pGMf 


~  (r2pV  )  +  r2  pV  V  -  -Li 

at  <t>  ar  r  4>  p 


B 

r  JL±  -  pv  V 
P  r  $ 


ST  r2  r-2- T +  l/2p(V  L  +  V  l)  +  _i  (B  Z  +  B  z) 
ot  y  -  1  r  <t>  2[i  r  <t> 


2  X  W  1  /q  2,-2) 


+  JZ-  r2v  — +  1/2  p(V  +  V  2) 
or  ry-1  r  4> 


ST(tV  +  3T  -  %  -  V,» 


ST  <*Br>  -  » 


(la) 

(lb) 

(lc) 

(ld) 

(le) 

(lf) 


The  various  constants  appearing  in  the  above  equations  are:  the 

rate  of  specific  heat  y,  the  magnetic  permeability  p,  the  gravitational 
constant  G,  and  the  solar  mass  Ms< 

Because  of  a  great  difference  of  the  two  methods,  however,  It  is  not 
possible  to  have  the  same  Initial  and  boundary  conditions  for  this 
comparison.  In  our  MHD  simulation  a  radial  velocity 

Vr  -  300.0  +  313.4  sin2  (Rt)  (2) 
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is  assumed  and  maintained  at  r  =■  0.5  au  for  all  t>0,  which  provides  a 
valid  solution  beyond  the  critical  point.  The  constant  Q  Is  chosen  so 
that  2n/0  »  25.4  days,  the  solar  rotation  period.  A  discussion  of  the 
finite  difference  method  of  solving  the  MHD  equations,  as  well  as  the 

initial  state  of  the  solar  wind,  may  be  found  In  Stelnolfsen,  et  al. 

(1975).  Computations  are  carried  out  to  10  AU  by  Dryer  and  Smith 
(private  communication).  We  ran  the  computation  until  the  first  shock 
advanced  well  beyond  a  distance  of  10  au.  Figure  5  shows  the  V-R 

relationship  at  T  =  2161.89  hours  (~  90  days). 

The  HA  method  requires  the  following  parameters:  A,  B,  i ^ ,  T,,  VQ, 
a  *  (R(a)  -  R(b))/(R(a)  -  R(c));  details  are  given  In  Hakamada  and 
Akasofu  (1982).  Figure  5  includes  the  V-R  relationship  and  the  R- 1 

relationship  for  the  following  set  of  parameters  at  T  =  90  days  =  2160 
hr: 

A  *  0.15  Vo  “  **80  km/sec 

B  -  0.85  a  »  0.7  exp  {-2  •  10"5  (jAtj-423)2} 

**  30000  hr  for  0  <  |At|  <  423  hrs 

m  550  hr  a  *■  0.7  for  423  hr  <  At 

It  is  quite  likely  that  a  number  of  sets  of  the  parameters  can  make  a 
reasonable  fit.  However,  the  non-uniqueness  Is  not  the  problem,  since 
what  matters  Is  the  V-R  and  R-t  curve  which  can  reproduce  accurately  the 
MHD  solutions.  In  spite  of  the  fact  that  an  exact  comparison  Is  not 
possible  (because  of  differences  of  the  origin  and  of  a  minor  difference 
of  the  time  variations  of  solar  wind  speed  at  the  source) ,  the  agreement 
between  the  two  curves  is  reasonable. 
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Figure  5.  Comparison  of  the  velocity  V  -  distance  R  relationships  obtained 
by  the  HA  method  and  a  one-D  MHD  simulation.  It  shows  also  the  density  and 
the  magnetic  field  magnitude  as  a  function  of  R  obtained  by  the  HA  method;  the 
results  from  the  one-D  solution  are  essentially  identical  to  those  obtained  by 
the  HA  method. 


An  important  point  to  make  here  is  that  it  is  relatively  easy  to 
modify  the  R-t  curve  in  such  a  way  that  the  resulting  V-r  curve  resembles 
quantitatively  a  solution  of  the  MHD  equations.  Taking  into  account 
various  uncertainties  involved  in  obtaining  the  MHD  solutions,  the 
agreement  between  the  two  V-R  curves  should  be  considered  satisfactory  at 
the  present  time.  If  Improved  one-dimensional  solutions  or  two-,  three- 
dimensional  solutions  become  available,  the  modification  could  be 
similarly  made.  If  space  probe  observations  could  eventually  provide  an 
accurate  V-R  relationship,  our  V-R  curves  would  also  be  improved  by  the 
same  procedure. 

In  concluding,  it  is  interesting  to  consider  why  a  crude  kinematic 
method  can  reproduce  reasonably  well  the  MHD  solution.  This  is  because 
the  solar  wind  flow  is  a  relatively  simple  gas  flow,  compared  with,  say, 
motions  of  the  terrestrial  atmosphere.  In  spite  of  the  fact  that  the 
collisionless  solar  wind  plasma  has  a  variety  of  complicated  internal 
processes,  there  is  no  major  internal  heat  source  and  boundaries  which 
alter  radically  the  initial  flow  speed,  except  in  the  vicinity  of  the 
planets  which  have,  however,  trivial  obstacles  in  the  heliospheric 
scale.  As  a  result,  the  R-t  relationship  is  well  confined  between  the 
two  curves  which  connect  the  maxima  and  minima  in  the  original  R-t 
curve.  That  is  to  say,  the  solar  wind  flow  has  a  significant  kinematic 
tendency.  This  tendency  has  recently  been  well  illustrated  by  Burlaga  et 
al.  (1983)  who  showed  that  the  solar  wind  speed  observed  at  Voyager  1  (~ 
8  au)  can  reasonably  be  inferred  from  Helios  (~  0.5  au)  observation, 
except  at  the  time  of  the  passage  of  the  shock  waves.  At  the  present 
time,  we  have  been  finding  even  a  better  method  to  improve  the  agreement. 


* 


A  Dynamo  Theory  of  Solar  Flares 

Most  of  the  modern  flare  theories  have  been  developed  on  the  basis 
that  the  magnetic  energy  stored  in  the  solar  atmosphere  is  suddenly 
released  by  magnetic  field  reconnection  or  by  current  Interruption. 
However,  there  has  been  little  discussion  of  how  the  stored  magnetic 
energy  is  generated.  Regardless  of  any  specific  models  of  solar  flares, 
it  is  important  to  treat  the  flare  phenomenon  In  the  context  of  an  energy 
generation-dissipation  system  as  a  whole.  On  the  basis  of  this  thought, 
we  have  proposed  a  dynamo  theory  of  solar  flares.  The  reason  for 
emphasizing  the  dynamo  process  is  two-fold.  The  first  is  to  propose  that 
temporal  variations  of  energy  dissipated  during  a  flare  follow  closely 
those  of  the  dynamo  power  which  can  supply  the  required  power  of  10  erg 

s-t  and  total  energy  of  10  erg.  This  is  in  contrast  to  the  prevailing 

view  that  the  energy  for  solar  flares  should  be  stored  in  the  coronal 
magnetic  field  prior  to  flare  onset.  The  second  is  to  formulate  the 
electrodynamic  coupling  of  the  photosphere-chromosphere-corona  system, 
powered  by  the  dynamo  process,  in  which  solar  flares  are  a  consequence  of 
enhanced  dynamo  power  output.  Thus,  we  emphasize  both  energy  generation 
and  energy  dissipation  associated  with  solar  flares. 

The  idea  of  powering  solar  flares  by  the  dynamo  process  in  the 

partially  ionized  chromosphere-photosphere  region  has  been  discussed 
earlier  by  several  workers,  including  Sen  and  White  (1972),  Heyvaerts 
(1974),  Obayashi  (1975),  Akasofu  (1979).  Certain  aspects  of  these  models 
are  briefly  reviewed  here.  Sen  and  White  (1972)  suggested  that  the 

dynamo- induced  Hall  current  leads  to  Joule  heating  in  the  photosphere  and 
that  the  electric  field  is  anti-parallel  to  the  convection  velocity  in 
the  photospheric  dynamo  region.  They  suggested  that  a  two-stream 
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instability  (Farley,  1963)  can  occur  in  the  Hall  current  to  trigger  a 
flare  in  the  photosphere.  Hayvaerts  (1974)  presented  a  photospheric 
dynamo  model  in  which  the  coronal  currents  are  entirely  field-aligned  and 
are  driven  by  horizontal  photospheric  motions.  Following  Alfven  and 
Cariqvlst  (1967),  Heyvaerts  (1974)  assumed  that  the  coronal  field-aligned 
currents  can  be  Interrupted  by  an  instability  (unspecified)  and  thereby 
release  the  stored  magnetic  energy  to  produce  flares.  Obayashi  (1975) 
proposed  a  flare  model  similar  to  Hayvaerts  (1974)  except  the  suggestion 
that  the  release  of  the  stored  coronal  energy  is  Initiated  by  a  collapse 
in  the  magnetic  field  and  the  formation  of  an  X- type  reconnection  line  in 
the  corona.  Recently,  Akasofu  (1979)  emphasized  that  the  dynamo  process 
in  the  photosphere  can  supply  the  required  power  (1029  erg  s-^)  and  the 
total  energy  (10  erg)  and  suggested  that  disruption  of  coronal  field- 
aligned  currents  driven  by  the  photospheric  dynamo  can  lead  to  an 
electric  potential  structure  which  accelerates  the  current-carrying 
electrons  to  produce  flares. 

The  proposed  solar  flare  process  to  be  developed  in  the  present 
paper  has  evolved  from  the  above-mentioned  dynamo  flare  models  (Ren  and 
White,  1972;  Heyvaerts,  1974;  Obayashi,  1975;  Akasofu,  1979).  The 
important  differences  between  the  proposed  model  and  those  earlier  models 
lie  in  (i)  solar  flares  are  treated  self-consistently  and  quantitatively 
as  a  consequence  of  the  electrodynamic  coupling  of  the  photosphere- 
chromosphere-corona  system  and  (ii)  the  onset  of  flares  is  quantitatively 
Identified  as  a  result  of  the  development  of  field-aligned  potential 
drops  due  to  the  loss-cone  constriction  effect  (Knight,  1973)  on  enhanced 
upward  field-aligned  currents,  Independent  of  any  Instabilities. 


Consider  now  the  photosphere-chroraosphere-corona  coupling  system 
which  Is  magnetically  connected  by  arch-like  (two-dimensional)  field 
lines  anchored  in  the  photosphere,  as  schematically  Illustrated  in  Figure 
6.  A  horizontal  motion  of  the  neutral  gas  in  the  photosphere  and  the 
lower  chromosphere  across  the  feet  of  the  arch-like  field  lines 
constitutes  a  dynamo  process.  For  simplicity,  the  dynamo  region  is 
assumed  to  consist  of  a  shear  flow  of  horizontal  neutral  wind  Vn  as 
indicated  in  Figure  6.  The  corona  is  located  in  the  upper  part  of  the 
arch  and  is  primarily  a  reactive  load  (capacitive  and  inductive).  The 
conjugate  chromosphere-photosphere  region  is  primarily  a  resistive  load 
(hereafter  referred  to  as  the  conjugate  region). 

As  the  suggested  horizontal  wind  blows  in  the  photosphere  and  the 
lower  chromosphere  (hereafter  referred  to  as  the  dynamo  region),  dynamo 
currents  are  generated  and  flow  through  the  whole  system  (the 
photosphere-chromosphere-corona),  both  along  and  across  magnetic  field 
lines.  The  power  generated  by  the  dynamo  consists  of  power  dissipated  in 
the  dynamo  region  and  power  output  to  the  load.  The  output  power 
delivered  to  the  load  is  partially  dissipated  in  the  conjugate  region; 
the  remaining  power  is  in  part  consumed  to  accelerate  current-carrying 
particles  along  magnetic  field  lines  and  in  part  stored  in  the  corona. 

Suppose  that  the  velocity  of  the  wind  in  the  dynamo  region  increases 
from  a  very  low  value.  As  the  dynamo  power  increases,  the  Intensity  of 
the  field-aligned  currents  will  be  enhanced.  However,  the  upward  field- 
aligned  current  density  is  limited  by  the  ’loss-cone'  effect  on  the 
precipitating  coronal  electrons,  due  to  the  converging  magnetic  field 
which  exerts  upward  mirror  forces  on  the  precipitating  electrons.  The 
'loss-cone'  is  defined  by  a  downward  facing  solid-angle  centered  along 
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Figure  6.  A  schematic  Illustration  of  the  photoaphere-chromosphere-corona 
coupling  system.  The  magnetic  field  B0  has  an  arch-llke  configuration  as 
Indicated  by  solid  curves.  The  dynamo  region  Is  located  around  the  boundary 
between  the  photosphere  and  the  lower  chromosphere.  Field-aligned  current  J( 
and  the  cross-field  current  are  Illustrated  by  dash-dot  curves.  VQ  Is  the 
neutral  wind  velocity.  E-  Is  the  dynamo  electric  field. 


the  magnetic  field.  A  particle  whose  pitch  angle  lies  Inside  the  loss- 
cone  will  reach  the  chromosphere- photosphere  level.  Once  they  reach  this 
level,  they  will  not  be  able  to  mirror  back  up  the  field  lines  due  to 

charge-neutral  collisions.  Only  those  electrons  Inside  the  loss-cone  can 
carry  the  upward  field-aligned  current  to  and  from  the  dynamo  region. 
Thus,  the  loss-cone  effect  Imposes  a  limit  on  the  upward  field-aligned 
current  density  (Knight,  1973;  Fridman  and  Lemalre,  1980).  To  exceed  the 
loss-cone  current  limit,  field-aligned  potential  drops  are  required  to 

accelerate  electrons  and  thereby  reduce  their  pitch  angles  Into  the  loss- 
cone.  The  required  field-aligned  potential  drops  can  be  distributed 
along  field  lines  and  supported  self-consistently  by  the  double  layer 

process  (for  reviews  see  FHlthammar,  1978;  Kan  and  Lee,  1981). 

As  the  accelerated  (upward  current-carrying)  electrons  descend  down 
along  the  legs  of  the  arch  toward  the  dynamo  region  and  its  conjugate 

region,  they  collide  with  the  gases  there,  Ionizing  and  exciting  them. 
The  resulting  optical  emissions  are  identified  with  the  optical 

manifestation  of  flares  during  the  rise  phase  excluding  the  Impulsive 

phase  (Kahler  et  al.,  1980). 

The  most  important  aspect  of  our  proposed  mechanism  is  that  the 

development  of  a  flare  follows  closely  the  temporal  variation  of  the 
power  output  from  the  dynamo  process.  Flare  onset  corresponds  to  the 

moment  when  the  upward  field-aligned  current  density  exceeds  the  loss- 
cone  limit  and  a  significant  potential  drop  Is  forced  into  existence 
along  the  magnetic  field  lines,  accelerating  the  current-carrying 
electrons  toward  the  lower  chromosphere.  The  maximum  epoch  of  a  flare 
corresponds  roughly  the  time  when  the  dynamo  power  output  maximizes.  Our 
theory  has  been  published: 


Kan,  J.  R. ,  S.-I.  Akasofu,  and  L.  C.  Lee,  A  dynamo  theory  of  solar 
flares.  Solar  Phys.,  84,  153,  1983. 

In  our  theory,  the  velocity  shear  Is  crucial  for  flare  production, 
and  It  Is  Important  to  verify  our  theory  on  the  basis  of  observations  of 
velocity  shear.  It  has  been  very  fortunate  that  Harvey  and  Harvey  (1980) 
have  already  examined  the  relationship  between  locations  of  the  velocity 
shear  In  the  photosphere  and  of  solar  flares.  Figure  7  shows  an  example 
of  their  observation  of  a  solar  flare  on  April  28,  1978. 

A  Simulation  Study  of  the  Formation  of  a  Bl-polar  Magnetic  Structure 

Emergence  of  magnetic  flux  through  the  photosphere  and  below  has 
been  discussed  extensively  In  terms  of  magnetic  buoyancy  (Parker,  1955, 
1979;  Tandberg-Hanssen,  1967;  Gilman,  1970;  SchUssler,  1977,  1979; 

Moffatt,  1978;  Meyer  et  al.,  1979;  Spruit,  1981).  The  dynamics  of 
magnetic  flux  has  mostly  been  discussed  In  terms  of  the  pressure  balance 
equation  pex  *  pjn  +  B£/8n,  where  pex  and  p^n  denote  the  plasma  pressure 
external  and  Internal  to  a  flux  tube,  respectively.  However,  such  an 
equation  may  be  applicable  only  for  a  situation  In  which  the  flux  tube 
rises  as  a  whole  along  an  extensive  horizontal  length  (much  longer  than 
the  extent  of  Its  vertical  motion).  It  will  be  hardly  applicable  for  a 
situation  envisaged  by  Babcock  (1961),  In  which  the  flux  tube  rises 
locally  with  a  Q-shaped  geometry  by  a  convective  motion,  such  a 
distortion  of  a  flux  tube  encounters  a  strong  Lorentz  force  which  tends 
to  restore  the  Initial  straight  geometry. 

Recently,  Galzauskas  et  al.  (1982)  and  Wallenhorst  and  Topka  (1982) 
showed  that  a  "complex  of  activity"  is  maintained  by  fresh  Injections  of 
flux  which  forms  and  disappears  In  rapid  succession  as  closely-spaced  but 
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Figure  7.  Boulder  region  1092  on  April  28,  1978,  Each  frame  covers  an  area 
of  256  by  256  arc  seconds  with  north  to  the  left  and  east  to  the  top.  (top 
left)  Velocity  image  with  black  representing  approach  and  white  recession, 
(top  right)  Velocity  shear  image  with  x's  indicating  flare  sites.  (bottom 
left)  Magnetic  field  with  white  representing  a  field  toward  the  observer  and 
black  away  from  the  observer.  (bottom  right)  Intensity  image. 


distinct  bl-polar  units;  most  of  the  magnetic  flux  disappears  in  less 
than  one  rotation  and  principally  in  situ,  without  evidence  of  signi¬ 
ficant  spreading.  If  magnetic  buoyancy  is  the  cause  for  the  formation  of 
a  Q-shaped  geometry  of  the  magnetic  flux  tube,  and  if  it  is  powerful 
enough  to  overcome  the  counteracting  Lorentz  force,  there  would  not  be 
any  other  force  to  submerge  it.  Therefore,  it  is  difficult  to  explain 
why  the  magnetic  flux  can  disappear  without  spreading. 

Therefore,  we  suggest  that  a  convective  motion  is  more  important 
than  magnetic  buoyancy  as  the  direct  cause  of  the  Q-shaped  magnetic  flux 
tube  thus  arisen  is  submerged  by  the  counteracting  Lorentz  force  as  the 
convective  motion  decays.  The  importance  of  the  convective  motion  was 
discussed  earlier  by  Weiss  (1966).  More  recently,  SchUssler  (1979) 
considered  vertical  motions  generated  by  magnetic  buoyancy.  Meyer  et  al. 
(1979)  examined  buoyant  magnetic  flux  tubes  in  supergranules.  Parker 
(1982)  examined  effects  of  flux  tubes  on  convection.  The  importance  of 
the  convective  motion  on  this  subject  was  also  recently  emphasized  by 
Raadu  (1983). 

On  the  basl3  of  the  above  consideration,  we  examined  the  formation 
and  disappearance  of  a  bi-polar  magnetic  structure  from  a  horizontal 
magnetic  field  structure  by  a  two-dimensional  MHD  simulation  method  by 
assuming  Initially  well-defined  convective  cells  which  are  allowed  to 
evolve  as  they  interact  with  the  magnetic  field. 

The  interaction  between  the  photospheric  gas  and  magnetic  flux  may 
be  described  by  the  standard  MHD  equations: 

^  x  (v  x  b)  +  \v2B  (1) 

o  t 
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(2) 


B)  x  B  +  g 

d  t  p  4np 

For  the  reasons  given  in  the  above,  magnetic  buoyancy  is  ignored  as  a 
first  approximation  by  assuming  that  the  mass  density  p  is  constant. 

We  choose  for  our  model  a  characteristic  length  L  **  50,000  km,  and  a 
characteristic  velocity  VQ  ■  0.1  km/sec  in  the  first  simulation  and  VQ  = 
1  km/sec  in  the  second  case.  Such  a  choice  of  the  values  of  L  and  V0, 
together  with  the  magnetic  diffusivlty  X.  “  10"^  km2/sec"^,  is  equivalent 

a  t 

to  choose  the  magnetic  Reynolds  numbers  of  the  order  of  X  x  10J  ~  x  10  , 
which  are  consistent  with  those  used  in  the  earlier  studies  (Weiss,  1966; 
Moffatt,  1978). 

The  MHD  equations  are  solved  by  an  explicit,  two-dimensional, 
formulation  of  the  Dufort-Frankel  method.  The  one-dimensional  Dufort- 
Frankel  method  is  given  In  Potter  (1973)  and  in  Richtmeyer  and  Morton 
(1967).  The  method  is  appropriate  for  parabolic  differential  equations 
and  is  numerically  stable  for  all  time  steps. 

First  of  all,  it  is  interesting  to  note  that  if  we  maintain  the 
velocity  pattern  within  the  convective  cell,  we  recover  the  earlier 
results  obtained  by  Weiss  (1966).  Our  result  is  shown  in  Figure  8  which 
may  be  compared  with  his  Figure  4. 

Let  us  assume  that  the  initial  velocity  pattern  is  allowed  to 
change.  Figure  9  displays  the  magnetic  and  velocity  field  configurations 
at  several  stages.  The  time  step  for  this  experiment  is  At  -  3536  sec  - 
1  hr.  The  maximum  fluid  velocity  is  observed  to  decrease  steadily  with 
Increasing  time  and  the  magnetic  field  lines  are  distorted  into  inverted 
"U"  contours.  Maximum  distortion  occurs  within  approximately  at  32At  or 
1.3  days.  The  most  important  aspect  of  this  particular  simulation  is 
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Figure  9  (1)(2).  The  deformation  of  a  uniform  magnetic  field  caused  by  a  convective  motion  at  different 
epochs.  The  time  step  A  t  is  about  1  hr.  In  this  case,  the  convective  motion  is  given  at  T  ■  0  and  let  it 
Interact  with  the  magnetic  field.  It  can  be  seen  that  the  convective  motion  is  reversed  after  about  1.3 
days  by  the  restoring  Lorentz  force. 


that  the  convective  motion  has  completely  been  reversed  by  the  distorted 


magnetic  field  after  the  maximum  distortion.  At  time  t  »  56At  or  about 
2.3  days,  the  magnetic  field  has  nearly  returned  to  approximately  its 
original  position.  As  this  process  of  "winding"  and  "unwinding"  the 
field  lines  continue,  the  maximum  absolute  velocity  tends  to  decrease. 

Comparing  Figure  3  and  9,  it  is  interesting  to  speculate  that  the 
actual  situation  around  sunspots  is  between  them.  The  magnetic  flux  can 
maintain  the  Q-shaped  geometry,  the  counteracting  Lorentz  force  tends  to 
restore  the  initial  horizontal  configuration,  forcing  the  Q-shaped  flux 
to  submerge  below  the  photosphere.  This  tendency  may  explain  the  new 
observations  by  Gaizauskas  et  al.  (1982)  and  Wallenhorst  and  Topka  (1982) 
that  most  of  the  magnetic  flux  in  a  complex  of  activity  disappears 
rapidly  in  situ  without  evidence  of  significant  spreading.  Our  results 
are  reported  ins 

Akasofu,  S.-I.  and  G.  Gislason,  A  simulation  study  of  the  formation 
of  a  bi-polar  magnetic  structure,  Planet.  Space  Scl.  (in 
press),  1984. 

It  is  our  plan  to  extend  our  simulation  study  of  sunspot  formation 
by  including  magnetic  buoyancy  and  other  factors. 

A  Study  of  the  Relationship  Between  Sunspot  Fields  and  Transient  IMF 
Variations 

The  configuration  of  the  magnetic  field  in  the  disturbed  solar  wind 
has  long  been  an  important  topic  for  solar,  cosmic  ray  and  magnetospheric 
physicists  (cf.  Gold,  1962).  This  topic  has  particularly  become  very 
Important  in  magnetospheric  physics,  since  it  has  been  found  that  the 
direction  of  the  interplanetary  magnetic  field  (INF)  is  an  Important 
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factor  In  causing  geomagnetic  disturbances.  Although  both  the  solar  wind 
speed  V  and  the  magnitude  B  of  the  IMF  are  also  important  and  are 
expected  to  be  large  In  the  disturbed  solar  wind,  the  north-south 
component  of  the  IMF  (Bz)  (accurately  speaking.  In  the  solar- 
magne tospherlc  coordinate  system)  plays  a  crucial  role  In  determining  the 
amount  of  solar  wind  energy  to  be  transferred  to  the  magnetosphere  (cf. 
Arnoldy,  1971;  Russell  and  McPherron,  1973;  and  Akasofu,  1981  and  the 
references  therein).  Specifically,  the  amount  of  the  transferred  energy 
becomes  very  large  when  the  IMF  has  a  southward  component  (~  -20  nT  for 
at  least  ~  6  hrs),  causing  intense  magne tospherlc  and  geomagnetic 

disturbances.  On  the  other  hand,  when  the  IMF  has  a  large  northward 
component  (~  +10  nT  for  several  hours),  the  amount  of  the  transferred 
energy  becomes  very  small,  resulting  In  only  a  weak  disturbance. 

The  present  large-scale  IMF  models  do  not  provide  such  large  values 
of  the  Bz  component  In  the  vicinity  of  the  solar  ecliptic  plane  (Parker, 
1958).  Thus,  as  In  the  past,  most  workers  assume,  implicitly  or 
explicitly,  that  the  Bz  component  is  associated  with  the  so-called 
'magnetic  tongue'  (Gold,  1962),  'magnetic  clouds'  (Klein  and  Burlaga, 
1982)  and  the  wavy  solar  current  disk  (Smith,  1981;  Akasofu  1981). 

The  magnetic  tongue  model  has  partially  been  supported  by  Pudovkin 
et  al.,  (1976,  1977,  1979)  who  claimed  that  the  degree  of  magnetic 

disturbances  can  be  predicted  by  the  north- south  component  of  the 
photospheric  magnetic  field  at  the  site  of  a  solar  flare,  implying  that 
the  polarity  of  the  IMF  can  be  predicted  from  the  flare  field. 
Specifically,  they  showed  that  If  a  aolar  flare  occurs  in  the  region 
where  the  photospheric  field  has  a  southward  component,  the  resulting 
geomagnetic  disturbance  tends  to  be  Intense,  vice-versa .  Their  study 
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Implies  thus  that  the  disturbed  solar  wind  generated  by  a  flare  carries 


out  the  magnetic  field  from  the  flare  region,  retaining  the  same  polarity 
of  the  field  at  least  to  a  distance  of  1  au.  In  a  somewhat  different 
context,  Lundstedt  et  al.  (1931)  demonstrated  that  the  solar  wind  speed 
tends  to  be  high  when  the  flare  field  has  a  southward  component,  while  no 
high  speed  wind  occurs  when  the  flare  field  has  a  northward  component. 
In  addition,  a  bi-directional  streaming  of  solar  cosmic  rays  has  also 
been  considered  as  a  supporting  evidence  of  the  magnetic  tongue  model 
(Palmer  et  al.,  1978).  As  noted  by  Lundstedt  (1982),  however,  one  of  the 
difficulties  in  this  study  Is  how  to  determine  accurately  the  orientation 
of  the  magnetic  field  in  the  vicinity  of  solar  flares.  This  is  because 
the  line  of  the  polarity  reversal  (or  the  so-called  'neutral  line')  tends 
to  run  in  the  north-south  direction,  since  sunspots  in  a  pair  tends  to 
align  in  the  east-west  direction.  Further,  solar  flares  take  place  most 
often  in  a  complex  sunspot  group.  As  a  result,  flare  'ribbons'  have,  in 
general,  a  very  complicated  geometry,  curving  in  various  directions. 
Therefore,  it  is  a  difficult  task  to  determine  a  single  value  for  the 
orientation  (say,  reckoned  clockwise  from  the  northward  direction).  The 
orientation  of  the  magnetic  field,  which  is  most  often  directed  across 
the  flare  ribbons,  depends  thus  on  where  it  is  measured. 

Tanaka  (1979)  found  that  there  is  a  high  degree  of  the  occurrence  of 
intense  solar  flares  in  the  configuration  called  "delta  spots"  which  are 
often  formed  by  non-paxred  spots.  In  delta  spots,  the  pair  thus  formed 
can  often  be  oriented  in  the  north-south  direction  or  in  the  east-west 
direction  (but  opposite  to  Hale's  law). 

North-south  oriented  delta  spots  are  of  great  Interest  here,  since 
the  line  of  the  polarity  reversal  tends  to  align  in  the  east-west 
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direction  and  thus  the  direction  of  the  magnetic  field  In  the  vicinity  of 
solar  flares  can  be  Inferred  with  reasonable  accuracy.  Tang  (1982,  1983) 
has  already  made  an  extensive  study  of  delta  spots  and  made  their  list 
(1969-1979);  here,  the  list  is  supplemented  by  including  recent  data 
(1980-1982).  Further,  In  this  particular  study,  Tange  et  al.  (1984,  In 
preparation)  examine  the  IMF  Bz  component,  instead  of  geomagnetic  indices 
which  can  be  Influenced  by  the  solar  wind  speed  V  and  the  IMF  magnitude 
B,  as  well  as  by  the  Bz  component. 

VI.  Results 


Delta  Configuration  and  Corresponding  IMF 

The  McMath-Hulbert  region  298  developed  Into  a  typical  delta 

configuration  and  passed  the  central  meridian  on  September  22,  1979 

(Tang,  1982;  see  Table  1).  Figure  10  shows  the  Kitt  Peak  magnetogram, 
the  Boulder  Hoc  filtagram,  the  Stanford  raagnetogram  on  September  22,  1979 
(Solar-Geophysical  Data,  November  1979  (No.  423)).  In  the  Kitt  Peak 
magnetogram,  white  portions  on  the  disk  Indicate  that  the  magnetic  field 
has  a  component  toward  the  earth  (positive),  while  dark  portions  have  a 
component  away  from  the  earth  (negative).  In  the  Stanford  magnetogram, 
solid  and  dashed  lines  Indicate  the  positive  and  negative  fields, 
respectively.  A  large  number  of  subflares  and  flares  were  observed  in 
the  region  between  September  16  and  24  (some  before  and  after  this 

period).  As  can  be  seen  In  the  Kitt  Peak  and  Stanford  magnetograms,  It 
was  the  only  prominent  active  region  around  the  central  meridian.  The 
positive  polarity  region  was  located  directly  north  of  the  negative 

polarity  region.  Such  a  polarity  is  described  as  N/S  polarity  for 
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SEPTEMBER  22.  1979 


Figure  10.  The  Kltt  Peak  magnetograra,  the  Stanford  magnetogram,  the  Boulder 
Hx  f iltergram  on  September  22,  1979  (Solar  Geophysical  Data)  and  the  IMF  B 
component  between  September  19  and  29,  1979  (ISEE-3). 


convenience  In  describing  the  IMF  orientation.  The  opposite  polarity  is 
described  in  the  S/N  polarity.  The  region  produced  one  3B  flare  on 
September  19,  five  2B  flares  (September  16,  17,  18,  19)  and  a  large 
number  of  IF,  B  flares  as  well.  A  medium  magnetic  disturbance  began  on 
September  19  and  lasted  until  September  22.  Another  medium  disturbance 
began  on  September  24  and  lasted  until  November  l. 

Figure  10  shows  the  IMF  Bz  component  between  September  19  and  29. 
It  is  clear  that  there  was  no  obvious  indication  that  the  N/S  polarity  in 
the  delta  configuration  extended  to  the  earth's  distance  during  or  a  few 
days  after  its  passage  across  the  central  meridian. 

We  have  just  begun  this  study  and  must  examine  a  large  number  of 
cases  before  any  conclusion  could  be  drawn.  If  we  could  find  a  simple 
relationship  between  the  solar  magnetic  field  orientation  and  the  IMF  Bz 
component,  the  degree  of  accuracy  of  predicting  the  occurrence  and 
intensity  of  geomagnetic  storms  will  increase  considerably.  Our 

preliminary  results  are,  however,  not  necessarily  very  promising. 
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